Decoupling heavy particles simultaneously 
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The renormalization group is extended to cases where several heavy particles are decoupled at the same time. 
This involves large logarithms which are scale-invariant and so cannot be eliminated by a change of renormal- 
ization scheme. A set of scale-invariant running couplings, one for each heavy particle, is constructed without 
reference to intermediate thresholds. The entire heavy-quark correction to the axial charge of the weak neutral 
current is derived to next-to-leading order, and checked in leading order by evaluating diagrams explicitly. The 
mechanism for cancelling contributions from the top and bottom quarks in the equal-mass limit is surprisingly 
non-trivial. 



1. SEVERAL LARGE SCALES 

When the renormalization group (RG) is ap- 
plied to a problem involving two or more large 
masses or momenta, the standard procedure is to 
deal with each asymptotic scale separately, start- 
ing with the largest. However, there are many 
cases of practical interest where this sequential 
approach will not do. For example, the top quark, 
the W ± and Z° bosons, and presumably the 
Higgs boson have masses in the same region, so in 
situations where more than one of these particles 
sets the scale, or there are momenta of similar 
magnitude, it is not reasonable to suppose that 
one of the corresponding logarithms dominates 
the others. Various techniques for dealing with 
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this situation are being investigated |1I2I3I4I5] , 

Our observation |6I7I8| is that the RG specifies 
simultaneous dependence on asymptotic scales 

Mi > Mi > . . . ^ M n > fixed scale ju , (1) 

for limits in which their logarithms grow large 
together: 

ln(Mi//x) /ln(Mj//i) = fixed, M u Mj -> oo . (2) 

The key new feature is the appearance of large 
logarithms 

HM./Mj) = ln(Mi//x) - HMj/n) (3) 

which are scale-invariant, i.e. independent of the 
renormalization scale /x. This precludes any use 
of the standard method |9I10| where a running 
coupling evolves from threshold to threshold and 
where, at each new threshold Mj, large scale- 
dependent logarithms are suppressed by choosing 
a new scale /i ~ Mj. 



Figure 1. Neutrino scattering from nucleons at 
low momentum transfers q. 



Figure 2. The lowest-order diagrams (two of them 
when fermion arrows are restored) responsible for 
mass logarithms induced by the neutral current 
©. Gluon propagators are shown as wavy lines. 
Straight lines denote propagators for heavy (h) 
and light (£) quarks. 



Instead, for each heavy particle hi (or large mo- 
mentum Q), we construct a scale- invariant run- 
ning coupling 



othi = ai H (la.(M 1 /ti), . . . ,\n.(M n /n)) 



(4) 



without reference to intermediate thresholds. All 
asymptotic dependence on M\ , . . . M n is carried 
by {othi] * = L ■ ■ ■ , «■}• Coefficients in the Taylor 
series in {a/^ } are determined by matching with 
the result for the sequential limit 



M, 



oo for fixed Mj +1 , . . . M n , j = 1,2,... (5) 



These ideas are tested for the classic example 
[9111112) of heavy quarks h being decoupled from 
the weak neutral axial-vector current 



J 
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at low momentum transfer q (Fig. . In partic- 
ular, we check that our results for simultaneous 
t, b decoupling are consistent with the sequential 
and rrifj ~ m t limits. 

Results are presented for three-colour quantum 
chromodynamics (QCD) with modified minimally 
subtracted (MS) renormalization at fixed scale p,. 

2. DIAGRAMMATIC ANALYSIS 

In perturbation theory, mass logarithms are 
first produced by the two- loop graphs of Fig. [21 



When both t and b are heavy, the result is |§] 



A subscript F refers to the number of flavours 
(F = 6 in this case); thus a 6 = (?g/47r 2 is the 
MS renormalized strong coupling for 6-flavour 
QCD. Eq. Q follows directly from Adler's two- 
loop renormalization Z-factor for the gauge- 
invariant axial- vector current in quantum electro- 
dynamics (QED). 

At this two-loop level, the manner in which mt 
and to;, are taken to infinity is immaterial: the 
limits can be taken separately or together, and 
as, cx.5, and need not be distinguished. 

However, we must be more careful when sum- 
ming all leading-order (LO) logarithms 



r+s+l 



ln r (mt//Li) ln s (m&//x), integer r,s ^ . 



For example, if mt and mb tend to oo together, 
as cannot enter the analysis, and the coupling 
of the residual F = 4 theory must be held fixed 
(not ae). 

Two types of diagrams have to be taken into 
account. First, there are the diagrams in Fig. [3 
which produce LO logarithms controlled by the 
RG; these depend on the strong coupling a>e of 
the original 6-flavour theory. Then the diagrams 
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Figure 3. Graphs responsible for leading-order mass logarithms in ghost-free gauges. Each gluon prop- 
agator is dressed with one-loop gluon and quark bubbles (all six flavours q) plus counterterms; these 
insertions sum to a geometric series, as shown. Tadpoles and other one-point gluon insertions are renor- 
malized to zero. Not shown are heavy-quark corrections to the light-quark propagator which produce 
finite corrections to make the residual singlet axial- vector current scale-invariant, as in Eq. Ijl2|l . 
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Figure 4. Leading-order matching condition for the /-flavour coupling a/ = g'j/'iir produced by decou- 
pling 6 — / heavy quarks h from the six-flavour theory with coupling — g^/Ait. 



of Fig. 0] are needed to obtain the LO matching 
condition expressing a$ in terms of the coupling 
ay of the residual theory with / light flavours. 



2.1. Preliminary Calculation 

First we check the LO amplitude for the 
top quark t to decouple from ij^j^t — ^7^756. In 
this case, the top mass m t acts as a regulator for 
the singlet component 



7 X! 97/*75«, t = u,d, s, c, b 



q =e 



of the current ^7^756 in the residual 5-flavour the- 
ory, with a renormalizing Z-factor produced by 
the diagrams of Fig. [3] in LO approximation. All 
dependence on a§ has to be eliminated in favour 
of «5, the quantity held fixed as mt — ► 00. So we 
apply Fig. 0] with / = 5, 

«5 = a6/ {l - (ae/37r) \n(m t /p)} , (8) 



and invert to obtain the matching condition 

a 6 = a 5 / {l+ (a 5 /37r) \a(m t /p)} ■ (9) 

LO 

After some work on the diagrams of Fig. |3J we 
obtain an asymptotic formula 

r M s ~ -(67 M 75% " (6/23) 2 7 M75 /ln(m t //2) (10) 

in complete agreement with the result of standard 
RG analysis |9I12| . The number 23 comes from 
the gluon self-energy insertions, i.e. from the one- 
loop factor 33 — 2/ in the beta function 



= -£( 33 - 2 /) - i!^( 153 - i9 /) 

+ 0{x 4 ) 



(11) 



of the residual f = 5 theory. The bare vertex 
7/i75 corresponds to the gauge-invariant singlet 
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operator normalised scale-invariantly |12I6) . 



q=i 



1= all / light flavours , (12) 
where / = 5 in Eq. ((TUj) . and where |14I12| 

7/W = J/ + ^3 ( 177 - 2 /)/ + °(^) ( 13 ) 

is the Callan-Symanzik function for the anoma- 
lous dimension of the singlet current. 

2.2. A Paradox? 

Suppose that both t and b are decoupled from 
^1^15^ ~ ^7/x75°- The residual theory has / = 4 
flavours, and all diagrams in Fig. shrink to the 
bare coupling 7^75 , which corresponds to the / = 
4 version of the invariant current (1121) . 

This case is interesting because of the role 
played by the equal-mass limit rrif, = m t , for 
which the answer to any order must vanish. To 
cancel the t term in i|l(J|) at TOf, = mt, a five- 
flavour b contribution has been suggested |9I12| : 

pcanccl „ (6/23) 2 7m75 

x {{Hm./fi))- 1 - {Hmt/fi))- 1 } . (14) 

However sequential decoupling of t and then b 
does not reproduce this result. When b is decou- 
pled directly from (67^756)5 in Eq. IjlOfl. standard 
RG theory gives [S] 

rJS> ~ (36/23) 7m75 

x{(251n(m 6 /j0) _1 - (23 hi(m t //i)) _1 } . (15) 

A pure five-flavour result is not obtained because 
the residual theory has four flavours - hence the 
factor 25 (= 33 - If for / = 4) in the b term. 

In fact, there is no paradox here: the limits 
mt ~ m;, and m t nib are extreme cases of the 
general limit To understand what is happen- 
ing, it is necessary to keep track of dependence 
on the scale-invariant logarithm ln(m t /mf,). 

2.3. Simultaneous t, b Decoupling 

The leading logarithms of diagrams in Fig. [3] do 
not depend on external momenta or light-quark 



masses me, so we set them equal to zero. The am- 
plitude can be written as a Euclidean integral over 
the momentum Q carried by each dressed gluon 
propagator, with self-energy insertions II(— Q 2 )- 



r t,b «6 

r " 5 To 2^ 7 " 75 



ds 



{1 




s(l 



dQ 2 



. Q 2 {!- n 
s){Q/m h fy X 



>)}' 



1 

h=b ' 



(16) 



The gluon integral has been restricted to the rel- 
evant asymptotic region \Q\ > m^. It can be fur- 
ther restricted by noting that, as t and 6 decou- 
ple, the factor in square brackets tends to 1 for 



mi 



< 



\Q\ ^n%t and is negligible otherwise: 



F " 5 l=o 2^ 7 " 75 



dQ 2 



Q 2 (i-n(-g 2 ))^ 



(17) 



In this Q-range, the gluon self-energy amplitude 
is 



n(-Q 2 ) 



23a 6 . m Q 2 _ 



a 6 ^ m^ 



(18) 



lo 127T 6ir y z 

so we are left with a simple logarithmic integral: 



r t,b 6a 6 
^ To 23^ 75 



mi , m 
23 In — - 2 In 



(19) 

fj, fl 

The last step is to match ae directly to 0:4 
(Fig. H for h = t,b) 

a 6 == a 4 / [l + (a4/37r){ln(m b //2) + ln(m t / p,)}] , 

(20) 

and eliminate a$ from l|19|) . This yields the com- 
plete set of LO logarithms 

6a 4 

A* 



r 



— + 2 In — 

fl (J, J 



(21) 



and hence the general asymptotic result [H] 



r t.b 36 
l ils ~ 23 W5 



{231n^ + 21n^V 
l u u J 



(22) 
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This answer replaces (|14f) and covers all cases: 

1. It vanishes for m t = m^. 

2. It reproduces the correct sequential result 
l|15|) when t is decoupled and then b. 

3. It works for intermediate cases as well: 
1 < \n(m t / ft) / ln(rrib / ft) < oo 

3. GENERALIZED RG 

Beyond LO, it is very difficult to identify all 
contributions "by hand" , so the RG becomes es- 
sential. Our aim |6I7I8| is to extend the RG to 
take all large scale-invariant logarithms into ac- 
count, and so generate asymptotic expansions for 
the simultaneous limit J5J: LO, NLO (next to 
leading order), NNLO (next to NLO), and so 
on. Since no choice of scale can eliminate all 
large logarithms, we simply fix the MS scale p, 
as {rrih} — > oo, as Witten ^H] did for a single 
heavy quark. 

Witten has all asymptotic dependence on rrih 
carried by a RG invariant running coupling 5/,. 
For (say) t decoupled from six- flavour QCD, a t is 
defined by the equation 

ra t 

ln(m t /p)= dx{l-S 6 (x)}/(3 6 {x), (23) 

J a e 

where Sf denotes the flavour Callan-Symanzik 
function for mass renormalization: 

8 F (x) = -2x/n + 0(x 2 ). (24) 

Then a§ is eliminated by a matching condition 
|16j relating it to the residual coupling a$ which 
is held fixed as mt — > oo; (similarly, light-quark 
masses have to be matched). This sequence can 
be continued by decoupling the b quark 

ln(m 65 //i)= f b dx{l-S 5 {x)}/f3 5 (x) (25) 

Ja 5 

with «5 matched to on, and so on. 

We begin by reformulating the matching pro- 
cedure to make RG invariance manifest. Once 



again, let t be the heavy quark, and let mt be its 
scale-invariant mass |15) : 

ln(m t /fi) = / dx/p 6 (x) . (26) 

J OS 

Then RG invariance implies the existence of a 
matching function J- defined by 7 

hx(fh t /p) = f dx/fi b {x)+T^ 5 (a t ). (27) 

The result of eliminating at from Eqs. I|26|l and 
H27fl is the all-orders matching condition relating 
ae to as . Similarly, there is a matching function Q 
for light-quark masses. For sequential decoupling, 
there corresponds a sequence of matching func- 
tions (^6^5, <?6^5), (-^5^4, £?5-mi) and so on. 
The T and Q functions are NNLO and NNNLO 
respectively. 

Now consider a simultaneous limit, e.g. the de- 
coupling of t, b. We seek scale-invariant running 
couplings a t and aj, which are defined without ref- 
erence to quantities a§ and m(, 5 associated with a 
b threshold, but which tend to Witten's couplings 
at and 5f, in the sequential limit J5|. 

Our prescription [S] is to start with Witten's 
formula Q23JI for the heaviest quark, and then use 
scale-invariant logarithms to generate the others. 
For t, b decoupling, we define at and at as follows, 

\n(m t /p) = / dx{l - 5 6 (x)} / (3 6 (x) , 

J a e 

\n(m t /m b ) = f *dx{l - 5 5 (x)} / /3 5 (x) , (28) 

J at 

where both m t and m^ are renormalized masses in 
the original sia;-flavour theory. The integrands of 
Eq. I|28|l are chosen such that, in the sequential 
limit, a t and a b agree with a t and a\> to NLO 
accuracy. 

The definitions of scale-invariant mass and 
matching functions also have to be generalised. 
We define invariant masses nit and rh~b for the 
simultaneous limit by the formulas 

ra t 

\n(m t /p) = / dx/f3 6 (x) , 

J a e 
rat 

\n(mt/m b ) = / dx//3 5 (x), (29) 

J a b 
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and then construct matching functions J- and Q 
for 6 — > 4 flavours, e.g. 



\n{m b /n) 



dx/P 4 (x) +T 6 ^4{a t ,a b ) . (30) 



The matching condition between ae and <x\ in the 
simultaneous t, b limit is found by eliminating at 
and a b from Eqs. and lj3"U|) . 

All {m/j} dependence in the leading asymptotic 
power of an amplitude is carried by the running 
couplings a t , a bl .... Thus, a power series expan- 
sion in a f , a b , ... yields coefficient amplitudes of 
the residual theory: 



Amplitude = 



^a^a 

ki... 



b ■ ■ 



k.C. 



This expansion remains valid for the sequential 
limit with the same coefficient amplitudes (to 
NLO — beyond that, there are calculable cor- 
rections due to T and Q functions). Hence these 
coefficients can be found by matching against the 
result for the sequential limit. 

For t, b decoupling from the weak neutral cur- 
rent, we found the result |6I8| 



X 1 + 



NLO 237T 

125663 6167 



(31) 



82800tt 3312tt 
Notice that Eq. I|28|) implies a t = a b for m t = m bl 
so both LO and NLO terms vanish in this case, 
as they should. Our definitions imply 

6 , ,36 



23tt 



\ctb 



at) 



23 



(32) 



-{ 2 31n^+21n^r 

so the LO term in (|31|) correctly reproduces the 
diagrammatic result (|22fl . 

The complete NLO result for the simultaneous 
decoupling of t, b, c from J^ eak is given at the end 
of 0. In that case, the definition 



In (rn b /i 



dx{l-S 4 (x)}/Pi(x) 



(33) 



is added to Eq. and Eqs. (O and (O are 

similarly amended so that ces is matched directly 

to Q!3. 
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